p propagation and dilepton production at finite pion density and temperature 
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We study the propagation properties of the p vector in a dense and hot pion medium. We 
introduce a finite value of the chemical potential associated to a conserved pion number and argue 
that such description is valid during the hadronic phase of a relativistic heavy-ion collision, between 
chemical and thermal freeze-out, where the strong interaction drives pion number to a fixed value. 
By invoking vector dominance and p saturation, we also study the finite pion density effects into 
the low mass dilepton production rate. We find that the distribution moderately widens and the 
position of the peak shifts toward larger values of the pair invariant mass, at the same time that 
the height of the peak decreases when the value of the chemical potential grows. We conclude by 
arguing that for the description of the dilepton spectra at ultra-relativistic energies, such as those 
of RH1C and LHC, the proper treatment of the large pion density might be a more important effect 
to consider than the influence of a finite baryon density. 
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I. INTRODUCTION 

One of the most salient features of the low-mass dilep- 
ton spectra in relativistic nucleus-nucleus collisions, from 
BEVALAC/SIS to SPS energies, is the enhancement in 
the production yields for invariant masses between 0.2 
and 1 GeV, as compared to proton induced reactions . 
Since dilepton final states are mediated by electromag- 
netic currents and these in turn are connected to vector 
mesons, dilepton pairs represent a prime tool to study the 
evolution of the dense and hot hadronic region formed in 
this kind of collisions. For low invariant masses, the vec- 
tor mesons involved are the p, tu and <f>. Among these, 
p plays an special role given that its lifetime is smaller 
than the expected lifetime of the interacting region and 
thus is able to probe different stages during the collision 
of heavy systems. 

The favored explanations, able to account for a great 
deal of the features of the measured low-mass dilepton 
spectra can be divided in two categories: the dropping 
p mass and the melting of resonances scenarios Q . The 
first of these, connected to the Brwon-Rho scaling conjec- 
ture Q and the decrease of the quiral quark condensate 
with temperature and baryon density, states that the in- 
medium p mass will sweep the entire low invariant mass 
region as the system cools down from its initially hot and 
dense state toward freeze-out. The second scenario 
states that the in-medium spectral densities of the p and 
its chiral partner, the a\ (1260) become broad and struc- 
tureless, merging into a flat continuum as the system 
approaches chiral symmetry restoration. 

An important ingredient for the success of the two 
above mentioned scenarios is the presence of finite baryon 
density effects (see for example Ref. Q)- In the first 
case, baryons act as the source of strongly attractive 
scalar fields. In the second case, interactions of the vector 



mesons with baryonic resonances, characterized by large 
coupling constants, overwhelm the relative strength of 
the interactions of p's with pions, despite the fact that at 
SPS energies, the relative abundance of the latter is five 
times that of the former. 

Nevertheless, at RHIC and moreover at LHC energies, 
the central rapidity region is expected to become baryon 
free with the relative abundance of pions being larger 
than at SPS energies. Consequently, at ultra-relativistic 
energies, it is important to include in the calculation of 
the dilepton spectrum the proper treatment of the large 
pion density, particularly during the hadronic phase of 
the collision. 

Recall that strictly speaking, pion number is not a con- 
served quantity and that pion decay is driven by all the 
relevant interactions, namely, strong, weak and electro- 
magnetic. However, the characteristic time for electro- 
magnetic and weak pion number-changing reactions, is 
very large compared to the lifetime of the system created 
in relativistic heavy-ion collisions and therefore, these 
processes are of no relevance for the propagation prop- 
erties of pions within the lifetime of the collision. As 
for the case of strong processes, it is by now accepted 
that they drive pion number toward chemical freeze-out 
at a temperature considerably higher than the thermal 
freeze-out temperature and therefore, that from chemi- 
cal to thermal freeze-out, the pion system evolves with 
the pion abundance held fixed 0, 13 ■ Under these cir- 
cumstances, it is possible to ascribe to the pion density a 
chemical potential and consider the pion number as con- 
served 0, 0| . In this context, the role of a finite pion 
chemical potential into a hadronic equation of state has 
been recently investigated in Refs. [II]]. The effects of a 
finite isospin chemical potential on the pion mass have 
also been recently studied in Refs. |T^ . 

The description of hadronic degrees of freedom belongs 
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FIG. 1: Feynman diagrams representing the pion self-energy 
at one loop. The wavy lines represent the p whereas the 
solid lines represent the pion. In the approximation where 
m p S> m w ,T, only diagram (a) contributes. 



FIG. 2: Feynman diagrams representing the p self-energy at 
one loop. The wavy lines represent the p whereas the solid 
lines represent the pion. 

effective Lagrangian. The theory thus obtained closely 
resembles scalar electrodynamics. We introduce a chem- 
ical potential associated to a conserved number of pions 
and construct the modifications to the pion propagator 
and 7r-/? vertex at finite density. In Sec. IIIII we use this 
effective Lagrangian to compute the pion self-energy and 
in Sec. II VI the p self-energy and from it, the modifications 
to its mass, width and dispersion curve at finite density 
and temperature. In Sec. |3 we compute the dilepton 
rate assuming vector dominance. We finally summarize 
and conclude in Sec. 

II. EFFECTIVE LAGRANGIAN 



to the realm of nonperturbative phenomena and there- 
fore has necessarily to rely on effective approaches that 
implement the dynamical symmetries of QCD. In a se- 
ries of recent papers |13L llij it has been shown that the 
linear sigma model can be used as one of such effective 
approaches to describe the pion propagation properties 
within a pion medium at energies, temperatures and den- 
sities small compared to the sigma mass. The sigma de- 
gree of freedom can be integrated out in a systematic 
expansion to obtain an effective theory of like-isospin pi- 
ons interacting among themselves through an effective 
quartic term with coupling a — 6(mJ/2/^), where 
and are the vacuum pion mass and decay constant, 
respectively. 

In this paper we extend the use of such effective de- 
scription to study the interaction of pions with the p 
vector, paying particular attention to the effects that a 
finite pion density introduce on the propagation proper- 
ties of p at finite temperature. We find the finite density 
and temperature modifications to the p mass, width and 
dispersion relation. By invoking vector dominance, we 
also study the effects that these modifications introduce 
on the production of e + e~ pairs near the p peak. 

The work is organized as follows: In Sec. [HJ we in- 
troduce the p in the description by gauging the original 



In Refs. [13| it has been shown that starting from the 
linear sigma model Lagrangian, including only meson de- 
grees of freedom and working in the kincmatical regime 
where the pion momentum, mass and temperature are 
small compared to the sigma mass, the two-loop pion self 
energy can be formally obtained by means of the effective 
Lagrangian 

C= 1 -(d^f- 1 -mW-^\ (1) 

where a = 6(to^/2/^) and the factor 6 comes from con- 
sidering the interaction of like-isospin pions in the vertex 

fff = -2i fe^ (5 ij S kl + 6 ik 5 jl + 6 tl P k ) . (2) 

In essence, the theory thus constructed and summarized 
by the effective Lagrangian in Eq. Q can be thought of 
as a theory for the effective coupling a that encodes the 
dynamics of low energy pion interactions. It can also be 
checked that Eq. reproduces the leading order modi- 
fication to the pion mass at finite temperature obtained 
from chiral perturbation theory |15| . 

In order to introduce the p field, we gauge the the- 
ory 0, 0| described by the Lagrangian in Eq. Q re- 
placing the derivative by the covariant derivative 
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FIG. 3: Thermal pion mass as a function of (a) T for different 
values of p, ranging from p, = to p, = 130 MeV and as a 
function of (b) p for different values of T ranging from T — 50 
to T = 150 MeV. 



FIG. 4: Thermal p mass as a function of (a) T for different 
values of p ranging from p = to p = 130 MeV and as a 
function of (6) p for different values of T ranging from T = 50 
to T = 150 MeV. 



given by 

.D^ = (0" -ig(f)<i>, (3) 

where we have introduced the ir-p coupling constant g. 
Also, by introducing the mass term and kinetic energy 
for the p field, the Lagrangian in Eq. becomes 

+ \myp»-\p^ v ■ (4) 

Let us now pause briefly to describe the formalism that 
allows us to introduce a finite chemical potential asso- 
ciated to a conserved pion number. To this end, let us 
further modify the Lagrangian in Eq. 10} , writing it in 
terms of a complex scalar field and regarding <j) and <fi* 
as independent fields 

+ Imyp^-^uP^ ■ (5) 

The effective Lagrangian in Eq. (JSJ resembles that of 
scalar electrodynamics with the photon field replaced by 
the massive p field. Invariance under global phase trans- 
formations 

0^0' = e~ lX (j) (6) 



leads to the conserved current 

J" = i(cj)*d^(p - <t> d"^*) - 2gp"(P*(f> , (7) 

and to the conserved charge N, that can be identified 
with the particle number, given by 

N = i J d 3 x(cj)*d (t) - (j) d (f>* + 2igp°<j>*<l>) . (8) 

Since N is a conserved charge, there exists a chemical 
potential fi conjugate to N so that the grand partition 
function is 

Z{fi,n) = Tr e-PiH-^N), ( 9 ) 

where j3 = 1/T is the inverse of the temperature T. From 
now on, let us work explicitly in the imaginary-time for- 
malism of thermal field theory. It is straightforward to 
write a path integral representation in Euclidean space 
of the grand partition function Z(/3,p). Care has to be 
taken by going through the Hamiltonian form since N 
depends on the time derivative of <fi and <fi* . After inte- 
gration over the conjugate fields it* = dC" /d(d°(f>*) and 
7r = dC" /d(d°(f)) one can check that in the exponent, 
there appear the combinations 

-^(^- 2 ^ + /x2 - m ')^' (10) 
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FIG. 5: Thermal p half width as a function of (a) T for dif- 
ferent values of p ranging from p, = to p, = 130 MeV and 
as a function of (b) p, for different values of T ranging from 
T = 50 to T = 150 MeV. 
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FIG. 6: p dispersion relation for (a) longitudinal and (6) trans- 
verse modes for p = 100 MeV for different values of T ranging 
from T = 50 to T = 150 MeV. 

The final outcome is that the introduction of a finite 
chemical potential translates into the substitution 

iu) n — > iuj n + p (15) 

both when this frequency appears in internal pion lines 
either in the Matsubara propagator or in the the ir-p 
vertex. Notice that these substitutions agree with the 
weel known result for the periodicity of the Matsubara 
propagator in the mixed representation given by 

A(/3-T,E;/i) = A(T, £;-/*). (16) 
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where r is the Euclidean time. Going to frequency space 
the Euclidean time derivatives get replaced by 



dr 

d(f>* 
dr 



(12) 



where the periodicity of the fields in the interval < 
t < (3 makes the (Matsubara) frequencies uj„ and oj n i be 
discrete and integer multiples of 27rT. The combination 
in Eq. I|1U[1 translates into a modification of the Matsub- 
ara pion propagator which now reads as (hereafter capital 
letters are used to denote four- vectors whereas lower case 
letters are used do denote the components) 



A(iuj n ,p; fx) 



1 



-(iuj n + A 1 ) 2 + P 2 + m 1 



(13) 



whereas the combination in Eq. <|ll|l goes into a modifi- 
cation of the n-p vertex which now becomes 

^7rp(P^,P'^^) = -ig{[~(iuJn + M),p] 

+ K&v + ^p']}, (14) 



III. 7T SELF-ENERGY 

The effective Lagrangian in Eq. 10) describes the in- 
teractions between pions and p's as well as among pi- 
ons themselves. The one-loop diagrams for the pion and 
p self-energies are depicted in Figs. 1 and 2. For the 
pion self-energy, the diagrams contain terms with inter- 
nal p propagators. Notice that since m p m x , T, these 
diagrams are strongly suppressed at finite temperature 
compared to those made out exclusively of internal pion 
lines and can thus be neglected. This approximation 
goes along the lines of the reasoning used in Refs. 01 
where the heavy internal sigma modes were systemati- 
cally pinched to construct the effective vertices and prop- 
agators leading to the effective Lagrangian in Eq. QJ. In 
this scheme, the pion and p self-energies decouple and 
the former gets modified only by the diagram in Fig. la 
leading to a momentum independent correction of the 
pion mass. We can thus take one step further going be- 
yond naive perturbation theory and adopt a resummation 
scheme for the pion self-energy to look for the modifica- 
tion of the pion mass beyond leading order. This can be 
implemented by writing the pion self-energy explicitly as 
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FIG. 7: Dilepton production rates for (a) longitudinal and 
(b) transverse modes for T = 150 MeV, k = 50 MeV and 
p = 0, 100 MeV. 



Equation (|17f) represents a self consistent relation for 
the temperature and density dependent quantity IIo. 
This is the well known resummation for the superdaisy 
diagrams which constitute the dominant contribution in 
the large- TV expansion |l8( of the Lagrangian in Eq. Q . 
The solution to Eq. (|17|) is given by the transcendental 
equation 
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n=l 



i\Jm\ + n \ cosh(n/i/T) 



T 



(18) 



Figure 3 shows the behavior of the pion thermal mass 
77T.Tr = W m 2 + Fin. as a function of (a) T for different 
values of p and (b) as a function of p for different values 
of T. We use the values = 135 MeV, f n = 93 MeV. 
From Fig. 3, we notice that tti,,- grows monotonically with 
both T and p. 



the vacuum pieces can be reabsorbed into the redefini- 
tion of the bare masses and couplings For what 
follows, we will concentrate on the matter contributions. 

For a massive vector field, the tensor structure of its 
self-energy can be written in terms of the longitudinal 
P£ v and transverse P^ v projection tensors 



W u = F(K)P£" + G(K)P% V . 



(21) 



In Minkowski space, Pj, L are given by 



puu _ 

pij _ 
T — 
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L — 



p° i = p4? = 



5 ij - k l k 3 /k 2 
-9 



■ K''I\";K- 



T ■ 



(22) 



From the relation between the self-energy, and the full 
D 1 *" and bare Dq U p propagators, we have, also in 
Minkowski space 



iD" v 



L 



T 



K 2 



F K 2 



G 



~m^K 2 



(23) 



In order to identify the coefficients F and G of the lon- 
gitudinal and transverse projectors, we take k along the 
z-axis. Thus, in Minkowski space, their expressions are 



IV. p SELF-ENERGY 

The explicit expression for the one-loop p self-energy 
depicted in Fig. 2 is given by 

™, = _ 2 TV f d3p {^-K^P v -Kn 
9 ^ J (2tt) 3 (P 2 + ml)[{K ~ P) 2 + ml] 

^ J (2/t) 3 P 2 + ml ' 



(19) 



where, according to the discussion in Sec. [HI the internal 
pion momentum and the p momentum are 

P" = [-(iu; n + p),p} 

K» = (-iw.k). (20) 
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(24) 



where II 03 and II 11 are obtained from Eq. (|19fl with the 
analytical continuation 



fc 



(25) 



that give the retarded functions and that can be per- 
formed after carrying out the sum over the Matsubara 
frequencies. 

The sum over frequencies can be obtained by standard 
techniques 0, |2(j ■ Considering the effects of a finite 
chemical potential, the expressions of interest are 
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T^A(iw n ,p)A(i(o;-u; n ),|k-p|) = 

n 

TyjcjnA ( to n , p) A (i (uj - |k - p|) = 



s1 S ± 4 ^-pi 



is 2 [1 + n(s 1 E p -fi) + n(s 2 -B|k- P | + /«)] 



E 

SlS 2 =± 



sis 2 (si^ p -//)[! + n{siE p - p) + n(s 2j E| k _ p | + p)] 



4£^E| k _ p | 



(iw - sx-Ep - S2^|k- P |) 



(26) 



where 



n(x) 



3 P\x\ _ 



(27) 



is the Bose-Einstein distribution, £" p 



-^■Ik-pl = ^/(k — p) 2 + m 2 and the function A is defined 
in Eq. ifBjl . 

Using Eq. 126(1 into Eq. 1(19(1 and by means of the an- 
alytical continuation in Eq. 1(25(1 . the real and imaginary 



= yjp 2 + rh 2 , parts of II 03 and II 11 are given by 
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Ren 
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8^1 dp vk) 



2(fc 2 + fc 2 ) 
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9 2 fl 



32tt VA: 3 

-fcp+afc 



dE p [4p 2 fc 2 - (A; 2 - fc 2 + 2/c £p) 2 ][l + 2n(£ p - k - fi)}6(k ) 



dE p [4p 2 k 2 - (fc 2 - fc 2 - 2fc £: p ) 2 ][l + 2n(£ p + fc + p)]6{-k )\ 



, . 4fcp + [fc 2 ± 2fc £ P ] In 

Z/C J) V 

(fc 2 - fc 2 ± 2fc £ p ) 2 - 4fcV 
4fc 3 p 



fc 2 - fc 2 ± 2fc -E p - 2fcp 



fc 2 - fc 2 ± 2k E p + 2kp 



In 



fc 2 - fc 2 ± 2fco^p - 2fcp 



fc 2 - k 2 ± 2fc £ p + 2fcp 



(28) 



where the function a is defined as 
/ 4rr? 2 

a= f-w^)- (29) 

It is easy to check that in the limit p — > 0, Eqs. 1(281) re- 
duce to the corresponding expressions found in Ref. [16| • 
Figure 4 shows the behavior of the p thermal mass 
obtained as the solution for fco of either 

fc o-^-^{S(fc^=!i)}- ' (30) 

as a function of (a) T for different values of p and (b) as 
a function of for different values of T. We use the value 



<? 2 /47r = 2.93 as determined by the p width in vacuum. 
For fc = there is no distinction between transverse and 
longitudinal modes and thus both Eqs. ((30(1 lead to the 
same solution. For every value of (T, /i), the solution 
is computed with the corresponding value for to w (T, p). 
From Fig. 4, we see that the thermal p mass grows mono- 
tonically with both T and p: and that the growth is larger 
for larger values of p. 

Figure 5 shows the behavior of the p decay rate (or 
half width) obtained from either 

/ F(m p (T, t x),k = 0) \ 
im \G(m p (T,p),k = 0) [ 

2m p (T,p) 
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as a function of (a) T for different values of p and (b) as 
a function of p for different values of T . Again, for k = 
there is no distinction between transverse and longitu- 
dinal modes. This can be shown analytically from the 
explicit expressions of ImF and ImG in Eqs. (|28|l which 
for k = yield 

ImF(fc o ,fc = 0) = ImG(fc , k = Q) 

~ tew -?r 

x [1 + 2n(fc / 2 + M)] ■ (32) 

For p — 0, Eq. l|32l) coincides with the corresponding 
expression in Ref. [lq . 

From Fig. 5a we see that for a given value of p the p 
width increases monotonically with temperature. From 
Fig. 5b we observe that the width reaches a minimum at 
a finite value of p and that this value increases as the 
temperature increases. This behavior can be understood 
if we recall that when the density increases so does the 
pion mass and thus the phase space available for the de- 
cay products of p narrows. The situation changes when 
the density is large enough so that the increase in the 
mass of the p becomes steeper than the growth in the 
pion mass, widening the phase space available for the 
decay process. 

Figure 6 shows the dispersion relation for (a) longitu- 
dinal and (b) transverse p modes for different values of 
T and p = 100 MeV. The main difference between the 
curves in each set can be attributed to the increase of the 
p mass with both p and T. 



V. DILEPTON RATE 

The electromagnetic current can be identified with the 
underlying quark structure of hadrons. For invariant 
masses below the charm threshold, this current can be 
decomposed as 

2 1 - 1 

J M em = 3«7^ - - -S7 M s . (33) 

Using the SU(3) quark content of hadrons, this current, 
being vectorial in nature, can be in turn identified with 
the current constructed out of the vector mesons p, u 
and 4>. 

3?* = (34) 

This is the well known conjecture named vector domi- 
nance model (VDM) . Since the pion electromagnetic form 
factor is is almost totally dominated by the p for invari- 
ant masses below 1 GeV, a simplified picture to study 
the decay of this electromagnetic current into low mass 
dilepton pairs is to consider that the current in Eq. I|34|) 
is totally dominated by the p. A further simplification 
stems from considering the spectral density of p as a sim- 
ple pole located at its peak which in turn dictates that 
the coupling of p to the electromagnetic current is em 2 p / g. 

Under the assumption of VDM and p saturation, the 
expressions for the thermal dilepton rates from longitudi- 
nal and transverse p modes are given, respectively, by fl6| 



E+E_ 



dR 



d 3 p+d 3 p- 
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d 3 p+d 3 p- 




(q • k) 2 



k 2 



-ImF 



(M 2 -m 2 - RcF) 2 + ImF 2 - 1 



2M Z 



(q-k) 2 



-ImG 



(M 2 



ReG) 2 + ImG 2 - 1 



,(35) 



where f>£ = (E+,p+) and Pi = (£L,p_) are the 
positron and electron four momenta, respectively, = 
P+ + Pi, <2 M = P+ - Pi, A/ 2 = k 2 - k 2 is the pair 
invariant mass, ujl.t are the longitudinal and transverse 
p modes dispersion relations and we have neglected the 
electron mass. 

Figures 7 and 8 show the dilepton production rates as 
functions of the pair invariant mass M for q ■ k = for 
fixed values of T = 150 MeV and k for two values of 
p = 0, 100 MeV for (a) longitudinal and (b) transverse 
modes. Figure 7 considers an small value of k — 50 MeV 
and Fig. 8 a larger value k = 250 MeV. In both cases we 
can see that the effect of the finite chemical potential is 
to moderately widen the distribution and to (also mod- 



erately) displace the position of the peak toward larger 
values of M. The most significant effect is however the 
lowering of the peak for finite p which is in agreement 
with the analysis of Ref. [2l| . 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have considered the effects of a finite 
pion density on the propagation properties of p mesons 
at finite temperature. The p has been introduced by 
gauging an effective Lagrangian obtained from the lin- 
ear sigma model in the kinematical regime where the 
pion mass and temperature are small compared to the 
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sigma mass. The finite density is described in terms 
of a finite pion chemical potential associated to a con- 
served pion number. We have argued that such descrip- 
tion is important for ultra relativistic heavy-ion collisions 
at RHIC and LCH energies where the central rapidity 
region is expected to become baryon free. In this situa- 
tions, we expect that the influence of the dropping p mass 
or melting or resonances scenarios to describe the dilep- 
ton spectra, which rely on the effects of a finite baryon 
density, become less important than the effects of the 
expected large pion density. This is so particularly dur- 
ing the hadronic phase of the reaction, between chemical 
and thermal freeze-out when the strong interaction drives 
pion-number to a fixed value. 

We have found that the p thermal mass increases 
monotonically with both the temperature and the pion 
density. However, the p width as a function of p and for 
fixed T starts off by decreasing, reaching a minimum at 
a finite value of p. This behavior can be understood by 
noticing that as the density increases, the pion mass does 
too and the phase space available for the decay process 
p — > 7T7T narrows up to a -temperature dependent- value 
of p. From this value, the increase in the thermal p mass 
is stronger than the increase in the thermal pion mass 
and this effect produces a widening of the phase space 
for the process, thus producing the p width to rise. 

Under the assumption of VDM and p saturation, we 
have also computed the dilepton production rate as a 
function of the pair invariant mass. We have found that 
the finite pion density produces a moderate broadening 
of the distribution and a moderate increase of the posi- 
tion of the peak. The finite pion density also produces 
a decrease of the distribution at the position of the peak 
compared to the p — case. 

The problem posed by the renormalization of the- 
ories involving resummation is by no means a simple 
one. In fact, in recent years this subject has been 
much actively pursued in the literature (see for exam- 
ple Refs. p3. 12a. l24j h The consensus is that when the 
theory is vacuum renormalizable, it will still be renor- 
malizable after resummation. The solution, which can be 
formulated using different languages, has been shown to 
require that the counterterms needed for vacuum renor- 
malization in ordinary perturbation theory need also to 
receive the benefits of resummation and be considered 
self-consistently in the resummation procedure. 

For the purposes of our work, let us stress that the re- 
summation we have implemented for the pion self-energy 



corresponds to the tadpole approximation (see Sec III of 
Ref. |22|). In fact, the gap equation that renders the 
finite temperature and density corrections to the pion 
mass in our work, Eq. (TSJ, is identical to the thermal 
part of Eq. (24) of the above reference; we have just 
gone one step further, expanding the integrand as a se- 
ries which allows us to analytically integrate each term. 
In Eq. (24) of Ref. |22j, the renormalization of the vac- 
uum self energy and four-point function have been car- 
ried out on the physical mass-shell condition and have 
taken into account the self-consistency. In the language 
of this reference, this is so because in the tadpole approx- 
imation, the renormalized four-point function is constant 
and given just in terms of the original coupling constant 
on the mass-shell condition. 

Had we gone beyond the tadpole approximation, the 
situation would have not be that straightforward. In fact, 
it has also been shown in Ref. 22j that the self-consistent 
resummation scheme requires a more complicated renor- 
malization of the four-point function to be taken into 
account alongside the renormalization procedure for the 
self-energy. The effects of this certainly more complete 
analysis might be interesting with regard to further ther- 
mal modifications to the pion mass in situations where 
the original coupling constant of the <f> A theory is much 
larger that one but is certainly beyond the scope of the 
present work where our coupling constant a, as deter- 
mined by the scale of the interactions, set by the vacuum 
pion mass, is of order 1 (in fact, what matters, as shown 
in Eq. (|18|l is the effective coupling constant given by 
a/(47r 2 ) which is of order 0.1). 

Finally, in order to predict the final dilepton spectra at 
RHIC and LHC energies, the results found in this work 
have to be placed into a model for the evolution of the 
collision and also possibly to include the effects of other 
vector or axial vector mesons. All this is for future. 
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